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Abstract 

A first-principles-based model is developed to investigate the influence of lead vacancies on 
the properties of relaxor ferroelectric Pb(Sc 1 / 2 Nb 1 / 2 )03 (PSN). Lead vacancies generate large, 
inhomogeneous, electric fields that reduce barriers between energy minima for different polarization 
directions. This naturally explains why relaxors with significant lead vacancy concentrations have 
broadened dielectric peaks at lower temperatures, and why lead vacancies smear properties in the 
neighborhood of the ferroelectric transition in PSN. We also reconsider the conventional wisdom 
that lead vacancies reduce the magnitude of dielectric response. 



I. INTRODUCTION 



Relaxor ferroelectrics are intriguing materials that characteristically exhibit a broad 
frequency-dependent dielectric response [1]. Although relaxors were discovered more than 
fifty years ago Q] the precise origin of their peculiar dielectricity remains controversial. Var- 
ious models have been proposed (e.g., Refs. DDDDDDQ and references therein) all 
of which include quenched, randomly distributed, internal electric fields and/or interacting 
polar microregions. 



Experimentally, the introduction of vacancies into relaxors induc es g reater broadening of 
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12j|. Other Pb- vacancy 



the dielectric response; i.e. it enhances relaxor character 
induced effects include: shifting of the dielectric peak to lower temperature; a decrease in 
the magnitude of the dielectric response; and a more "diffuse" phase transition 
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The mechanisms responsible for vacancy-induced characteristics are mostly unknown, 
but understanding them would illuminate the fundamentals of relaxor behavior, especially 
because defects are always present in real samples. One reason for this paucity of knowledge 
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is that first-principles studies of (point) defects in ferroelectrics are scarce [13, 

In this manuscript, we investigate the effects of Pb- vacancies on the properties of a relaxor 
ferroelectric from first principles. Our model system is Pb(Sci/2Nbi/2)03 (PSN) rather than 
the prototype relaxor Pb(Mg!/ 3 Nb2/3)03 (PMN), because the former is easier to model - 
i.e. finite-temperature first-principles-based models already exist for PSN j^, 0| but not for 
PMN. Note that, unlike PMN which retains relaxor character down to OK (i.e. exhibits no 
ferroelectric phase transition), chemically disordered PSN exhibits relaxor behavior between 
~ 400K and ~ 368K, then transforms to a rhombohedral ferroelectric phase We 
focus on static properties, because studying the dynamics of relaxors from first-principles 
is not currently feasible. Our study clarifies the origin(s) of the vacancy-induced effects 
mentioned above, and suggests that the common wisdom that Pb-vacancies reduce the 
magnitude of dielectric response should be reconsidered. 

This manuscript is organized as follows: Section II describes the method. Section III 
presents and discusses results. Section IV is a brief summary. 
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II. METHOD 



Our method is a generalization of the first-principles-derived effective Hamiltonian of 
Ref, HHQto invest.ga.te p—es with che m ical disorder on their B-sites and 
Pb-vacancies on their A-sites, by writing the total energy as a sum of two terms: 



£PSN-DV({M, r] H , {r?/}, {<7j}, {s k }) 

= £ P sn-d({uJ,?7h, {?7/}, {<Jj}) + E^dui}, {s k }) , (1) 

where -Epsn-dv is the total energy of the Pbi-zd^Sci^Nbi/^Os alloy (PSN-DV), □ rep- 
resenting a vacant A-site; .Epsn-d is the total energy of chemically disordered PSN solid 
solutions (PSN-D); E yaic gathers the explicit energetics associated with Pb-vacancies; u, 
is the (B-centered) local polar distortion in unit cell i (proportional to the local dipole 
moment); t]h and {rjj} are the homogeneous and inhomogeneous strain tensors 2(|, respec- 
tively; CTj = —1 or +1 if there is a Sc or Nb cation, respectively, at the B-lattice site j of 
the Pbi_ x n a; (Sci/2Nbi/2)03 alloy Finally, {s^} characterizes the amount and distribution 
of A-site vacancies, i.e. Sk=0 or 1 if the fc'th A-site is occupied by a lead atom or vacant, 
respectively. 



For i? PSN _ D , we use the analytical expression proposed in Ref. |17|, which includes a 
term of the form — J2i ■ ei[{crj}] , where Z* is the Born effective charge associated with 
the local distortion, and e; is the internal field at cell i that is caused by the surrounding B- 
site cation distribution [21[. 

For E vac , we use the following pertubative expression (first-order in Uj): 

E vac ({ui}, {s k }) = S k,i s k gki ■ Ui (2) 

k i 

where the sum over k runs over A-sites, and the sum over i runs over B-sites. The S^, 
parameters quantify how vacancies perturb the local distortions, is the unit vector joining 
the A-site k to the B-site-center of Uj. We truncated Sk,i at the first-neighbor shell. All 



the parameters of Eqs. (1) and (2) are derived from first-principles calculations 22l l23l \24\ 
performed on relatively small supercells (i.e., up to 40 atoms), including one (to obtain Sk^i) 
that exhibits a lead vacancy and a background charge of —2. 
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We solve this effective Hamiltonian by the Monte Carlo (MC) technique |25[ to simu- 
late properties of two different systems: 1) disordered Pb(Sci/2Nbi/2)03 without any Pb- 
vacancies, PSN-D; 2) disordered Pb n qrD q ns fSci /?Nbi 72)03, PSN-DV, which has a realistic 
5% vacancy concentration on its A-sites (loL V A 12 x 12 x 12 supercell (8640 atoms, 48 
A lateral size) is used for both systems. The variables Oj and are randomly chosen, and 
remain constant during the simulation. Temperature, T, is decreased in small steps from 
high T. We found it unusually difficult to get good statistics and, therefore, performed 10 
MC runs of 4 million sweeps each, at every temperature. We collect the supercell-average 
u of the local distortion at each sweep and classify it as visiting one of the eight rhombo- 
hedral, twelve orthorhombic, or six tetragonal energy minima by computing its projection 
onto each of these twenty-six high-symmetry directions. The visiting rate of a particular 
minimum is defined as the ratio between the number of sweeps in which that minimum is 
visited and the total number of sweeps. The Xa/3 dielectric susceptibility tensor is computed 
as in Refs. j^E^ that is by using the polarization fluctuation formula: 

(N Z*) 2 . 

Xa/3 = , „ [< U a Up > - <U a ><Up >\ (3) 

V e KB-L 

where < u a up > denotes the statistical average of the product between the a and (3 com- 
ponents of the supercell average of the local mode vectors, and where < u a > (respectively, 
< up >) is the statistical average of the a- (respectively, (3-) component of the supercell 
average of the local mode vectors. N is the number of sites in the supercell while V is its 
volume, ks is the Boltzmann's constant and e D is the permittivity of the vacuum. Note that 



a detailed derivation of Eq.(3) is provided in Ref. 27] 



III. RESULTS 

The simulated temperature at which the dielectric response is maximized, T m , is about 

I I 

850K for PSN-D and 794K for PSN-DV |28J]. Thus, our model qualitatively reproduces 
the previous finding that Pb- vacancies reduce T m jlU [ill Q ■ Orrr quantitative prediction 



of a 6.5 % reduction in T m , from the introduction of 5% P 



>vacancies agrees rather well 
and is also consistent with 



with the corresponding experimental reduction of 8.5% ^ 

the ~ 2.6% reduction of T m associated with 1.7% Pb-vacancies Q|; establishing the ac- 
curacy of simulations based on Eqs. (1) and (2). Note that Eq.(2) can be rewritten as 
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E vac ({ui},{s k }) = ~T,iZ*Ui ■ e vaCii [{s k }} , with e vaCii = -^Ek% s k g ki ; that is, e vaCii 
has the dimensions and physical meaning of an internal electric field arising from the Pb- 
vacancies. Our first-principles-derived values for the fields around a Pb-vacancy are of the 
order of 1.4 xlO 9 V/m at each neighboring B-site i. These fields are oriented towards 
the vacancy because Sk,% is positive and e vaC) i depends on the directional g& vectors; e va c,i 
is therefore rather large and inhomogeneous, and thus significantely opposes homogeneous 
ferroelectric {uj} ordering, which explains the vacancy-induced reduction of T m . 

Figure la displays the (u x ,u y ,u z ) supercell averages of local distortions as functions of 
the reduced temperature T/T m ; [ (u x ,u y ,u z ) is directly proportional to the spontaneous 
polarization, where the x, y, and z axes are chosen along the pseudo-cubic [100], [010], and 
[001] directions, respectively]. Both PSN-D and PSN-DV are predicted to have a low-T 
rhombohedral ferroelectric phase in which u x ~ u y ~ u z , in agreement with experiment 
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111 ] . Therefore, only the average component, (u) = l/3(u x + u y + u z ), is plotted in 
Fig. la. Also shown in Fig. la is (\u\)r/\/3, where (\u\)r is the mean magnitude of the 
supercell average of local distortions calculated from sweeps in which the system visits 
rhombohedral minima. Figure lb shows x ~~ which corresponds to one third of the trace of 
the Xa/3 dielectric susceptibility tensor - versus T/T m . In other words, and according to Eq. 
(3) and the definition of x, Fig- lb displays: 

(N Z*) 2 1 r 



X 



Ve k B T 3 



< u x + Uy + u 2 z > -(< u x > 2 + < u y > 2 + < u z > 2 ) (4) 



Predictions for (u) and x are quite scattered around T m , which reflects our difficulties 
in obtaining good statistics for either PSN-D or PSN-DV, as expected for systems with 
complex energy landscapes jaj. We smoothed the results for (u) by fitting them to a 7th- 
order polynomial in the region around T m (lines in Fig. la). We denote the result by (w)fit 
and recompute the susceptibility as 



*\2 



x» = ^rr(5W«-< fi >St), (5) 

where Oq is the OK cubic lattice constant. We obtain Eq. (5) from Eq. (4), by approximating 
V as Na^, and the averages as follows: 1/3 < u 2 + u 2 + u 2 >~ 1/3(|m|)|j and l/3(< u x > 2 
+ < u y > 2 + < u z > 2 ) ~ (w)flt |22|- The "smoothed" Xfit's are depicted with lines in 
Fig. lb, and will be used in combination with (w)fit, to illuminate the effects of Pb- vacancies 
in relaxors. 
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Interestin gly, Fig , lb indicates that, contrary to common wisdom and previously reported 



observations 
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11 1 , our predicted maximum value of the dielectric response does not de- 



crease when Pb- vacancies are introduced: our simulations indicate that xat is slightly larger 
in PSN-DV than in PSN-D at T = T m , because T m is smaller and (u)flt is further away from 
(\ u \)r/V3 hi PSN-DV, see Fig. la and Eq. (5). [Our predicted maximum values for x&t are 
about 34,000 and 38,500 for PSN-D and PSN-DV, respectively, which is close to the mea- 
sured value ~ 40,000 of Ref. 1(| for PSN-D]. This result, which appears to be at odds with 
measurements, is explained by the experimental observation of Ref. [3^ that Pb-vacancies 
in Pb(Sci/2Nbi/2)03 favor the simultaneous formation of pyrochlore phases (Pb2.31Nb2O7.31 
and Pb 3 Nb 4 0i3). Thus, dielectric measurements on a PSN-DV sample at T = T m should 
yield a smaller response than those on PSN-D; because, while the majority of the PSN-DV 
sample is perovskite, with a dielectric response that is at least as large as that for PSN-D, 
a smaller fraction is pyrochlore, which has a much smaller dielectric response 3^| . 

Figure lb also shows that our simulations successfully reproduce two experimental obser- 
vations: (1) the broadened dielectric peaks in PSN-D and PSN-DV, which are characteristic of 
relaxors 1[, and according to Eq. (5) and Fig la, arise from the diffuse phase transition asso- 
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ciated with (u)^; (ii) Pb- vacancy induced dielectric broadening ^0, 
width of the x&t peak is about 0.19 in PSN-D versus 0.22 in PSN-DV (in reduced units) [c.f. 
the corresponding predicted value of 0.13 reported in Ref. |9j for the "normal" ferroelectric 
system Pb(Zro.6Tio.4)03]. Note that vacancy-induced dielectric broadening only occurs below 
T m . This is particularly clear in the Fig. lb inset, in which the x&t of PSN-DV is rescaled 
by a T- independent constant to match the x&t dielectric response of PSN-D at T = T m . 

To clarify how PSN-DV differs from PSN-D only at temperatures below T m , Table I lists 
the internal energies of rhombohedral and orthorhombic minima (calculated at 5K); plus 
the visiting rates of these minima in both PSN-D and PSN-DV at T/T m ~ 0.93. The energy 
landscape of PSN-D is rather anisotropic; e.g. the energy difference between the ground state 
and the highest-lying rhombohedral minimum is 0.68 meV. Reference jj| demonstrated that 
such local symmetry breaking is caused by the anisotropy of the large internal electric 
fields (which have an average magnitude ~ 2 xlO 9 V/m; determined numerically) that arise 
from the charge-difference between Sc 3+ and Nb 5+ ions. This anisotropy results from 
our nanometric periodic PSN-D supercell being too small to contain a perfectly random 
B-site conation 3. Re S n,t s in Table I als0 indicate that, at T,T m , 0,3, PSN-D 
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preferentially visits the rhombohedral [-111] and [111] minima, i.e. the minima that have 
lowest energies, and that are connected by one of the lowest-in-energy orthorhombic minima 
(specifically, [Oil]). Such fluctuations lead to a (w)fit that is smaller than {\u\)r/v3, and 
thus [Eq. (5)], to a large dielectric response relative to normal ferroelectrics (which are 
already stuck in a single minimum at T/T m ~ 0.93 j^). As indicated by Table I, the 
e vac,i electric field associated with Pb-vacancies (see rewriting of Eq (2) above) brings the 
orthorhombic and rhombohedral minima closer to each other; similar to the finding of Ref. 
Q] that €i in PSN-D reduces barriers between energy minima. Thus, PSN-DV fluctuates more 
and/or visits more minima than PSN-D at temperatures below T m . This implies a smaller 



(it) at (Fig. la) and there 
observed experimentally 



ore naturally explains two vacancy-induced effects that have been 
lillllQ: a broader PSN-DV dielectric response (Fig. lb and Eq. 



(5)); and a more diffuse PSN-DV phase transition (PSN-DV "only" begins to stick in a single 
rhombohedral minimum at T/T m ~ 0.85 as indicated by the equality (w)fit = {\u\)r/\/?> at 
and below this temperature, Fig la, whereas, freezing in of PSN-D occurs at T/T m ~ 0.90). 

Results in Table I also indicate that the most frequently visited rhombohedral minima 
in PSN-DV are not the lowest-energy minima; rather they are those with small-in-energy 
rhombohedral-to-orthorhombic barriers (e.g., the [1-1-1] and [-1-1-1] minima that are con- 
nected by the [0-1-1] minimum, over an energy barrier of about 2.5 meV). This strongly 
suggests that relaxors with Pb-vacancies can be thought of as non-ergodic systems 

IV. SUMMARY 

In summary, first-principles-based calculations predict that the introduction of Pb- 



re temperature at which the dielectric response 



vacancies in Pb(Sci/ 2 Nbi/2)0 3 : (i) reduces t 

peaks, T m , consistent with experiment [h], 11, 3b (") does not intrinsically decrease the 
magnitude of the dielectric peak, contrary to common wisdom jiol . llll . Q]; (iii) broadens 
the dielectric response and enhances the diffuse character of the ferroelectric transition, 

nnn 

consistent with measurements [lQj, |11|, |12j (but only at temperatures below T m ). Item (i) 
above is the direct consequence of the relatively large internal inhomogeneous electric fields 
generated by the Pb-vacancies. We suggest that our disagreement with experiments in item 
(ii) is caused by parasitic pyrochlore phases in the experimental samples. Item (iii) origi- 
nates from the decrease in energy barriers that are caused by Pb-vacancies, which promote 
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enhanced fluctuations between different energy- minima below T m . Interestingly, our results 
regarding the effect of defects on energy landscape may also be relevant for modeling (and 
understanding) dynamical properties of relaxors. 
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TABLE I: Total energy per five-atom cell E, at 5K, and visiting rate v r at T/T m ~ 0.93 associated 
with the different rhombohedral (111) and orthorhombic (110) minima in PSN-D and PSN-DV. 
The zero of energy corresponds to the ground-state. The uncertainty is lower than 0.01 meV and 
0.1% for E and v r , respectively. 



Minimum 


E (meV) in PSN-D 


E(meV) in PSN-DV 


v r in PSN-D (%) 


v r in PSN-DV (%) 


[-1-11] 


0.38 


0.41 


6.3 


0.0 


[-111] 


0.00 


0.23 


72.9 


0.0 


[111] 


0.07 


0.00 


20.0 


0.0 


[1-11] 


0.68 


0.58 


0.0 


0.0 


[1-1-1] 


0.00 


0.28 


0.0 


37.2 


[11-1] 


0.67 


0.33 


0.0 


17.3 


[-11-1] 


0.65 


0.38 


0.0 


5.8 


[-1-1-1] 


0.20 


0.32 


0.0 


37.5 


[-101] 


3.05 


2.90 


0.3 


0.0 


[011] 


3.05 


2.80 


0.5 


0.0 


[101] 


3.45 


2.85 


0.0 


0.0 


[1-10] 


3.42 


3.42 


0.0 


0.0 


[10-1] 


3.48 


2.93 


0.0 


0.7 


[01-1] 


3.85 


2.99 


0.0 


0.2 


[-10-1] 


3.60 


3.12 


0.0 


0.2 


[0-1-1] 


3.05 


2.85 


0.0 


1.1 


[0-11] 


3.71 


3.29 


0.0 


0.0 


[110] 


3.69 


3.34 


0.0 


0.0 


[-110] 


3.54 


3.31 


0.0 


0.0 


[-1-10] 


3.58 


3.53 


0.0 


0.0 
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FIG. 1: Properties of PSN-DV (open dots, dashed lines) and PSN-D (filled dots, solid lines) as 
functions of T/T m . Dots in panel (a) show the supercell-averaged mean component of the local 
distortions ((it)). Lines represent fits to a 7th order polynomial ((u)&t)- In panel (a), we also show 
(squares) the magnitudes of the local distortions, (\u\)r, divided by \/3- Panel (b) displays one 
third of the trace of the dielectric susceptibility tensor directly obtained from MC simulations. The 
lines represent the Xfit results obtained from Eq. (5). The inset of Panel (b) displays x&t results 
that were rescaled to match Xfit of PSN-D at T = T m . 
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